Non-linear charge reduction effect in strongly-coupled plasmas 
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Abstract 

The charge reduction effect, produced by the nonlinear Debye screening of high-Z charges oc- 
curing in strongly-coupled plasmas, is investigated. An analytic asymptotic expression is obtained 
for the charge reduction factor (/ c ) which determines the Debye-Hiickel potential generated by a 
charged test particle. Its relevant parametric dependencies are analyzed and shown to predict a 
strong charge reduction effect in strongly-coupled plasmas. 

PACS numbers: 51.50+v, 52.20-j, 52.27.Gr 
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A basic aspect of plasma physics is the so-called Debye shielding of the electrostatic po- 
tential. This consists in the property of plasmas (or electrolytes Q), either quasi- neutral or 
non-neutral, to shield the electrostatic field produced by charged particles (to be also denoted 
as test particles) immersed in the same system. This result has fundamental consequences 
on plasma phenomenology, since it actually limits the range of static Coulomb interactions 
inside the Debye sphere, i.e., at a distance p < \d from the test particle, Ad being the 
Debye length. As usual here A^ = A^J , where the sum is carried out on all plasma 



species and Xd s = y attZ^n ' ^ anc ^ being respectively the s— species temperature 
and number density (the latter defined in the absence of test particles). In fact, when both 
the particles and the plasma are assumed non-relativistic, small-amplitude, stationary (or 
slowly time- and space- varying) , electrostatic perturbations generated by isolated test parti- 
cles, result effectively shielded in the external domain, i.e., at distances larger than the Debye 
length Xd. The renewed interest in this problem is particularly related to dusty plasmas 
or colloidal suspensions 

as 

which are characterized by the presence of a large fraction of 
highly charged particles (grains), i.e., having an electric charge Z^e with \Z^\ ^> 1. It is well 
established that the phenomenon of Debye shielding of the electrostatic potential generated 
by a slowly moving or stationary charged test particle, when occurring in strongly coupled 
ionized gases, manifests peculiar properties. These are produced, as a consequence of the 
nonlinear plasma response occurring at distances smaller than the Debye length Xjj, when 
the electrostatic potential (<E>) results locally such that $ = |e$| /T Q ~ 1 (or even ^> 1), $ 
being the normalized electrostatic potential evaluated at a suitable characteristic distance 
p from the position of the test particle r(t) and T D = Tj. The distance p Q can be identified 
with the radius of the local plasma sheath around a test particle. For spherically-symmetric 
test particles such a domain can be identified with the spherical subset of the configuration 
space, of radius p Q = ( J (average interparticle distance), iV = N oi being the density 
of the ion species, in which each test particle can be considered as isolated. Ionized gases 
(or electrolytes) can be classified according to the characteristic dimensionless parameters 
%o = Po/^d = g 1 ^ 3 and T = f3/x Q (where (3 = A ^ e \ D is denoted as the dimensionless elec- 
tric charge of the test particle), g and T being respectively so-called plasma and Coulomb 
coupling parameters. In particular, for plasmas the high-density requirement x z Q = g 1 
is satisfied by assumption and can be interpreted as the condition that the total number of 
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mutually interacting particles in a Debye sphere (expressed by 1/g) results 3> 1. Instead, 
the ordering of parameter V remains in principle arbitrary. Thus, the orderings T 1 and 
T ~ 1 (or r ^> 1) correspond respectively to so-called weakly and strongly-coupled plas- 
mas. Dusty and colloidal plasmas may be characterized by up to 10 4 — 10 6 with typical 
grain size smaller than 10~ 4 cm and with plasma temperature and the ion density T ~ 1 
eV, N Q ~ 10 s — 10 10 cm -3 . In this case the Coulomb parameter for a negatively- charged 
grain in the presence of the plasma sheath produced by hydrogen ions can result as large 
as T = (10~ 2 — 10~ 3 )Zd, while the dimensionless radius of the ion plasma sheath can be 
estimated x a = 0.01 — 0.05. Therefore, dusty plasmas typically result strongly coupled if 
Zd > 10 3 . For these plasmas it is important to be able to evaluate the electrostatic potential 
generated by dust grains which is expected to be strongly influenced by Debye shielding. 
It is well-known that, in general, this phenomenon occurs provided suitable physical as- 
sumptions are introduced. In particular, the plasma must be assumed appropriately close to 
kinetic Maxwellian equilibrium, in which each particle species is described by a Maxwellian 
kinetic distribution function carrying finite fluid fields [defined respectively by the number 
density, temperature and flow velocity (N, T, V)] . In the absence of test particles these fluid 
fields must be assumed slowly varying in a suitable sense, or constant, both with respect 
to position (r) and time (t). In this regard it is important to remark that the appropriate 
treatment of the plasma sheath surrounding each test particle is essential also for the validity 
of the mathematical model for the Debye screening problem, i.e., for the existence of classical 
solutions of the Debye screening problem, which do not exist when letting x Q = As a 
consequence the electrostatic potential depends necessarily on the dimensionless parameter 
x Q . In particular, in the domain x > x Q the dimensionless potential & Xo (x) generated by 



a dust grain of charge Z^e < 0, fulfilling the boundary conditions 
lim a; ^ 00 $ Xo (x) = 0, can be proven to satisfy the integral Debye-Poisson equation 



— — and 
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px POO 

dx'x' 2 + / dx'x' 

x 



S(x',x ), (1) 



yielding for x = x Q the constraint equation 

POO 

$ = T- / dx'x'S(x',x ). (2) 



Here <3> = $ Xo (x ), while in the case of a negatively-charged test particles in a quasi-neutral 
three species (electron, ion and dust) plasma the source term S(x, x D ) reads S(x, x D ) = 0(x — 



x Q ) &exp{|Z;| ®x (x) \ - £ e exp j-$ x . o (x) j - £ d exp |- \Z d \ $ Xo ( x )\ , Q(x-x ) being the 

i A 2 

Heaviside function. Moreover, £ s is defined by the ratio £ s = j^jf-- In the case in which 
the contribution due to the dusty species is negligible, and moreover \ZJ\ = 1, Tj = T e , 
the previous expression reduces to the customary value S(x,x Q ) = 0(x — x D ) sinh Q Xo (x). 
Another useful representation of Eq.Q can also be obtained by means of the transformation 



$ Xo (x) = $> exp{y(x,x )} , (3) 

and there results, by consistency with the integral DP equation ((TJ), y(x ,x ) = and 
y'(x, x )\ x=Xo = —a, with a = T/x & . This delivers for y(x,x D ) the integral equation 

y(x,x Q ) = y (x,x ) + F(x,x ), (4) 
y (x,x ) = -a(x - x ), (5) 

where F(x, x Q ) is the solution of the integral equation 

X 

F(x, x ) = J dx' (x -x') is (x', x ) - y(x', x Q ) 2 - (^'j x o) | , (6) 

Xo 

with S(x', x ) to be expressed in terms of y(x f , x Q ) by means of Eq.®. On the other hand, for 
finite x Q , the Debye-Hiickel asymptotic approximation for the electrostatic potential can be 
proven to hold when measuring the potential at a position r sufficiently far from the position 
of the test particle r(t). In fact this involves imposing that the so-called weak-field condition 
^ Xo (x) 1 be locally satisfied. As a consequence, the normalized electrostatic potential 
generated by a spherically-symmetric point particle of charge Ze can be approximated by 

where Ax = x — x , q = q(x ,(3) is a suitable dimensionless effective electric charge and 
3?£a (x) is usually known as DH potential. The latter identifies the external asymptotic 
solution which holds in the subset of the external domain, outside the Debye sphere, in which 
the weak-field condition is satisfied. The remaining notation is standard, thus p = |r — r(£)| 
is the distance from the point charge and x = p/\ D & [x Q , oo] is the corresponding normalized 
distance. The weak-field condition is manifestly locally satisfied both for strongly and weakly 
coupled plasmas, provided there results x > 1, with x suitably large (i.e., x 1). However, 
in the case in which there results either j:„ <C 1 and/or T ^> 1, i.e., for strongly-coupled 



plasmas, the electrostatic potential may actually decay on a scale much shorter than the 
Debye length and precisely on a scale suitably close to the boundary of the plasma sheath, 
i.e., for x Q < x <C 1. This conclusion is consistent both with previous analytic estimates 
based on the exact solution of the ID Debye screening problem \^ and direct numerical 
simulations of the 3D Debye screening problem (f| Q, 0, l|. It implies, however, that the 
DH approximation typically holds also at distances suitably larger than the boundary of the 
plasma sheath x Q and which extend also to values comparable to the Debye length, i.e., for 
. Finally, as a further consequence, the effective electric charge carried by the 
DH potential q may appear significantly reduced with respect to the case of weakly coupled 
plasmas 0, 1^1 0, Q]- To display this effect, it is convenient to represent the normalized 
effective charge of the test particle in the form q = /3f c , f c being the charge reduction factor. 
By dimensional analysis, it follows that f c must be of the form f c = / c (x , f3, £ s ), i.e., f c is 
function, to be determined, of the only dimensionless parameters of the problem, namely x Q 
and T (or f3) and moreover the ratios £ s (with s = e,i,d). It is well-known that for weakly- 
coupled plasmas q can be approximated by the normalized electric charge of the isolated 
test particle, i.e., q = (3, which implies f c = 1. Instead, for strongly-coupled plasmas, the 
value of charge reduction factor, determined either via one- dimensional analytic estimates 
Q or by m ea„s of a variety of — —on .netoods (see, for example, RflflQ), 
has been found to be notably smaller than unity, thus suggesting the existence of a possible 
significant charge reduction effect. Indeed the investigation of the effective interactions 
characterizing high-Z grains in plasmas has attracted interest in recent years especially for 



their role in dusty plasmas 
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121 ] . However, for strongly-coupled plasmas the precise form 



of the function f c is still unknown and, in particular, an analytic estimate of the effective 
charge characterizing the DH potential in strongly-coupled plasmas is not yet available. 

Besides being of obvious interest as a still unsolved mathematical problem, it is especially 
important from the physical standpoint to estimate f c as a function of the parameters x Q 
and T for strongly coupled plasmas characterized either by a high density and/or by the 
presence of a large fraction of highly charged particles (grains), such as dusty plasmas and 
colloidal suspensions 2 L 3]. 

In a previous work the Debye screening problem (DSP) has been formulated for an 
electron-hydrogen plasma in order to analyze the asymptotic properties of its solutions near 
the plasma sheath (i.e., for x suitably close to x D ). For this purpose the case of a strongly- 
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coupled plasma was investigated in which the parameters x Q and T satisfy the asymptotic 
conditions, to be denoted as strong- coupling ordering, T ~ ^> 1 with x a ~ 0(S k ), 
k = 0, land 5 an infinitesimal. The similar case of three-species plasma can be obtained 
directly. As a result, invoking for $ Xo (x) the representation (jH)), together with Eqs.(jIJ),(|HJ) 
and imposing the appropriate boundary conditions for § Xo (x), it follows that the integral 
on the r.h.s. of Eq.(j2J) can be estimated asymptotically to yield an asymptotic equation for 
the initial condition <l> , i.e., 

$ (x , P) S 4t In ( [r - $ 1 \ . (8) 



In particular, for an electron- hydrogen plasma in which the contribution to Ad due to the 
dusty species is negligible, there results \Z,-\ = 1 and & = 1/2 and this reduces (jSJ) to the 
expression given in Ref. In addition, neglecting higher-order infinitesimals of order 0(5 r ) 
(r > being a suitable real number), the following inequality can immediately be proven 



f c <c^(x ,(3)[l + 0(5 r )}, (9) 
c {a \x ,(3) =x ^ (x ,(3), (10) 

f3) being an asymptotic upper bound for the charge reduction factor. This result, 
as indicated below, is relevant to investigate the charge reduction effect in strongly-coupled 
plasmas. 

In this paper we propose a more accurate asymptotic approximation for f c , potentially 
useful to analyze the charge reduction effect of high-Z test particles taking place in strongly- 
coupled plasmas. In particular, we intend to prove that, in validity of the previous strong- 
coupling ordering, the charge reduction factor can be approximated in the form: 

f c ( Xo ,P) - c^(x ,P) = Xce^fixJ/P, (11) 
where, x c is defined so that x c > x Q and is provided by the root of the equation 

y'(x c ,x ) = -1- -. (12) 
x c 

Here $i* n *' ) (x) denotes the leading-order contribution to the internal asymptotic solution of 
the Debye-Poisson equation (0), which can be obtained from the representation 
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and the integral equation © and holds in an appropriate neighborhood of the plasma sheath 
^J. There results 

(x)-exp{y (x,x )}, 



while y'(x,x ) reads 



x 



y'(x, x Q ) = —a — a Ax + 2a log h 

Xq 

X 

+Tsr- / dx'exp | |Zj| $ exp {— a Ax'} + a Ax' H — a 2 Ax' 2 — 2a 

$o J I 2 



(13) 



(14) 



x' log Ax' 

X 



where Ax' = x' — x Q . Eqs.(JTTJ) and (fT2j) uniquely follow by imposing suitable matching 
conditions between the internal and external asymptotic solutions, (x) and $z"^(x), 
precisely 



r? ^ 

dx *' 1 y 



x,. 



X r : 



(9 ^ 

_$( ea: *)(' 

dx x ° 1 



(15) 
(16) 



It is immediate to prove that Eq. (jl6|) is equivalent to the condition of stationarity for the 
function fc a \x, x a , ft) 

i) 

0. (17) 



~Q^f C \ X T X Olft) 



Therefore Eq. ()17|) also furnishes the same asymptotic approximation for the charge reduction 
factor given by Eq. ffTTj) . The accuracy of the asymptotic approximation , ft) and of 

the upper bound c^ a \x Q , ft) can be tested by comparing them with the estimate of the 
charge reduction factor obtained from direct numerical simulations of the 3D Debye screening 
problem. For this purpose the function 



g(x,x Q ,ft) = xe x x °$ Xo (x)/ft 



(18) 



has to be evaluated numerically by solving the Debye-Poisson equation. The charge reduc- 
tion factor is defined as the limit 



f c (x Q ,ft) = lim xe x x °$ Xo {x) / 'ft. 



(19) 



In practice, for numerical estimates it can be approximated by f c (x , ft) = g(xf,x ,ft) with 
Xf to be suitably defined. It is important to stress that, in the limit 5 — > 0, one expects 
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x c — > 0, xj — > and moreover 



limc (6) (x ,/5) = hmf c (x ,P). (20) 

5^0 <5^0 



In addition, it must result also 



limc (fe) (x ,/5) = \imc {a) (x ,p). (21) 

Therefore, for 5 small enough (i.e., for T ^> 1 and/or x Q <C 1), it it actually possible 
to approximate f c (x Q , 0) both in terms of , /3) and Numerical simulations 

indicate, however, that the asymptotic estimate c^ b \x ,(3) actually converges rapidly to 
f c (x , (3) and furnishes therefore a more accurate approximation than the upper bound 
(3). Nevertheless, particularly thanks to Eq.(|21|). c^ a \x , (3) can still be used for order- 
of-magnitude estimates. 

In figures^the result of two numerical simulations are reported, for the case of an electron- 
hydrogen plasma and for = 1/2, which correspond respectively to (x Q , (3) = (0.08, 2.5) and 
(x , (3) = (0.15, 10). In both cases the matching point x c is found to be close to the boundary 
of the plasma sheath, which implies that the DH potential actually approximately applies for 
x > x c . The upper bound estimate (x a , (3) can be used for order-of-magnitude estimates of 
the charge reduction effect. Finally, the parametric dependencies of (3) with respect 

to (3 and x Q are displayed for selected values of the parameters. Figures El demonstrate that 
for typical values of the Coulomb parameter V and the radius of the plasma sheath x Q the 
effective electric charge of the DH potentially appear dramatically reduced. 

The present results are relevant for the investigation of strongly-coupled plasmas, such 
as dusty plasmas, in which a significant fraction of high-Z charged particles is present. 
In fact, the charge reduction effect has fundamental consequences on the Debye screening 
phenomenon, since it actually limits the magnitude of the Coulomb interactions produced by 
these particles. The estimate of the charge reduction factor f c , here obtained, is a necessary 
prerequisite for the descr ipt i on of a variety of kinetic processes (for example, condensation 



nption 

HQ 



collisional processes |lJ,ll5j, etc.), as well as electrostatic or electromagnetic instabilities 



in dusty plasmas (see for example |16l Il7||). 
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FIG. 1: Comparison between the asymptotic approximation of the charge reduction factor 
c( a \x ,(3), the numerical estimate f c (x , (3) and the upper bound estimate c^ a \x ,(3). Case a) 
corresponds to x Q = 0.08 and = 2.5 for which there results x c = 0.098, (x , (3) = 0.211, 
f c (x ,P) = 0.246 and (x , /3) ^ 0.33. Case b) corresponds to x Q = 0.15 and p = 10 for which 
there results x c ^ 0.171, d b \x ,/3) = 0.106, f c (x ,f3) ^ 0.119 and cW(x ,/3) ^ 0.19. 
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FIG. 2: Plot a) shows c^ a \x , T) as a function of with x Q = 0.1, 0.06, 0.024, 0.048). Plot b) shows 
c<°) (x , r) as a function x Q for = 10, 20, 40, 60, 80, 100. 
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